Physical Chemistry Chemical Physics

The molten globule (MG) state is an intermediate in the unfolding pathway of proteins, typically triggered by denaturing agents such as urea, extreme pH, high pressure, or heat. The microscopic details of such states are far from understood. Here we study the MG states in protein Hen Egg-White Lysozyme (PDB ID: 1AKI) using microscopic constant pH molecular dynamics (CpHMD) simulations and experiments across a wide pH range. We observe that the titratable residues act as key drivers of conformational fluctuations, promoting the emergence of MG states at extreme pH. These states display partial unfolding, and small global structural changes (< 7% deviation). Hydration around the fluctuating acidic residues shows reduced water density and weakened hydrogen bonding at low pH. At high pH, hydration around acidic residues increases relative to pH = 7, whereas hydration around basic residues decreases. The translational and rotational dynamics of hydration water also exhibit pronounced pH dependence: the translational diffusion coefficient (Dtrans) increases linearly with decrease in pH in acidic medium and increases linearly with increasing pH in the basic regime. The rotational diffusion (Drot) shows similar dependencies on pH except a break at pH {approx} 4 corresponding to acidic residue pKa values. Our results may be useful to identify ligand binding of lysozyme in extreme pH conditions.

There has been a surge in antibiotic resistance in recent years, making traditional antibiotics less effective against key pathogens. RNA has recently emerged as a potential target for antibiotics due to its involvement in crucial microbial functions. It is possible to expand the range of therapeutic targets by using RNA-based therapies, but it remains necessary to improve the molecular-level understanding of interactions between RNA and known and potential binders. The SAM-I riboswitch, which controls the transcriptional termination of gene expression involved in sulfur metabolism in most bacteria, is an excellent ligand target. Thus, understanding its behavior with and without ligand complexes would be very helpful for drug design applications. In this manuscript, we studied the interactions between the SAM-I riboswitch and its natural ligand, SAM, which controls riboswitch function, and compared those interactions to its interactions with the very similar small molecular SAH, which does not control riboswitch function, and to its interactions with a potential binder JS4, identified via virtual screening. From our simulations, we gain a deeper understanding of small molecule interactions with the SAM-I riboswitch. The results reveal how differently the small molecules (SAM, SAH and JS4) bind to and potentially induce conformational changes in the riboswitch. Our findings offer valuable insight into the molecular mechanisms underlying riboswitch RNA-ligand interactions for the design of more effective RNA-targeting therapeutics.

The conformational dynamics of a model cationic protein in water and in the presence of anionic sodium dodecyl sulphate (SDS) and cationic cetyltrimethylamonium bromide (CTAB) surfactants at different concentrations were investigated using all-atom molecular dynamics simulations. Free-energy landscapes constructed along principal components reveal a compact, well-defined native basin at 25 {degrees}C in water, whereas elevated temperature (100 {degrees}C) induces a broadening of the conformational space and the emergence of multiple metastable states. The presence of surfactants further modulates this behavior in a concentration-dependent manner. Cluster population analysis shows that SDS promotes a highly heterogeneous ensemble characterized by reduced dominance of the native-like cluster, while CTAB partially protects the protein from thermal denaturation at higher concentrations. Radial distribution functions demonstrate strong accumulation of SDS headgroups around the protein and pronounced insertion of SDS alkyl tails into hydrophobic protein regions, indicating direct hydrophobic destabilization and micelle-assisted unfolding. In contrast, CTAB exhibits weaker headgroup association owing to electrostatic repulsion and reduced tail-hydrophobic contacts, suggesting a less disruptive interaction mechanism. At high concentration, CTAB aggregates provide a structured hydrophobic environment that stabilizes the folded state and suppresses denaturation. Together, these results provide a molecular-level picture of how surfactant chemistry and concentration govern the conformational stability of a cationic protein, highlighting the dominant role of hydrophobic interactions in surfactant-induced denaturation at high temperature. Graphical Abstract O_FIG O_LINKSMALLFIG WIDTH=200 HEIGHT=89 SRC="FIGDIR/small/717321v1_ufig1.gif" ALT="Figure 1"> View larger version (24K): org.highwire.dtl.DTLVardef@dcf96aorg.highwire.dtl.DTLVardef@17acdc7org.highwire.dtl.DTLVardef@15bdc2borg.highwire.dtl.DTLVardef@1d39f3c_HPS_FORMAT_FIGEXP M_FIG C_FIG

In mammalian cells, lipid monolayers support the integrity of lipid droplets (LDs), organelles that function as storage for neutral lipids. Liver-targeting illnesses such as liver cancer interrupt normal LD metabolism and prompt changes in the chemical content of these organelles, which can have effects on structural and organizational behavior of the lipids. In LDs, liver cancer induces concentric crystalline phases of cholesteryl esters (CEs) and triglycerides near the NL-monolayer interface, which become more pronounced as CE concentration increases. Yet, there is little known about how this phenomenon may link to persistence of undigested LDs in liver cancer patients. To shed light on this, all-atom molecular dynamics simulations were used to model LD micropipette aspiration experiments and gain insight into the effect of CE concentration on partitioning, structural, and mechanical properties of LDs. We successfully model micropipette aspiration by application of constant surface tension laterally, which stretches lipid bilayers and monolayers as the magnitude increased. The results show increased phospholipid packing due to insertion of CE fatty tails into the monolayer. Increasing CE concentration induces a non-linear change in surface packing defects on the LDs, notable rigidification, and stiffness. Taken together, these insights improve our understanding of the physical properties at the LD monolayer-core interface during liver cancer progression.

Mycobacteria are responsible for causing severe illnesses like tuberculosis and leprosy in humans. Studying the mycobacteria cell envelope presents a significant challenge due to its intricate lipid compositions and structural variations and also its harmful nature in a typical experiment setting. In this study, we use all-atom molecular dynamics simulation to study mycobacterial inner membranes (MIMs). By incorporating different types of phosphatidyl-myo-inositol-mannosides (PIMs) and their glycoconjugates such as lipomannans (LM) and lipoarabinomannans (LAM) lipoglycans, we have constructed both symmetric and asymmetric membrane systems to study the MIM structure and dynamics under varying compositions of each lipid type. Our results show that the phospholipid/PIM-rich inner leaflet remains a stable, fluid bilayer, and the outer leaflet structure and dynamics are heavily governed by lipoglycan surface density. Importantly, as LM/LAM concentration increases, the polysaccharide chains shift from flexible, membrane-lying orientations to a compact brush-like state aligned with the membrane normal. This crowding significantly reduces the solvent-accessible volume and limits direct interactions between LM/LAM sugars and the outer leaflet surface. Furthermore, we observe that high lipoglycan presence in the outer leaflet slows lipid diffusion across the entire bilayer, demonstrating a dynamic coupling between the two leaflets. By resolving these LM/LAM sugar-level dynamics and their impact on membrane-wide properties, this study provides a molecular framework for future MIM modeling and simulation with various (peripheral) membrane proteins to better understand how the MIM functions as a regulated physical barrier and a platform for mycobacterial virulence.

Intrinsic tryptophan fluorescence is widely used as a sensitive reporter of protein conformational dynamics, yet the molecular origin of its temperature-dependent modulation remains unclear. Here we investigate the conformational dynamics of Trp134 in bovine serum albumin (BSA) using molecular dynamics (MD) simulations, free-energy calculations based on umbrella sampling and WHAM, quantum mechanical (QM) calculations, and QM/MM approaches. MD simulations show that the global structure of BSA remains stable while temperature induces a gradual population shift from the Ia+ to the Ia- rotamer. The corresponding free-energy landscapes reveal that this shift arises from subtle changes in basin stability and transition barriers along the rotameric coordinate. In contrast, standalone QM calculations on isolated tryptophan predict different energetic trends, highlighting the sensitivity of rotamer stability to electronic-structure treatments and environmental effects. QM/MM calculations partially reconcile these differences by incorporating the protein environment. Together, these results suggest that temperature reshapes the rotamer free-energy landscape of Trp134, leading to population shifts that modulate intrinsic tryptophan fluorescence in proteins.

Formation of the {beta}-amyloid (A{beta}) plaques is a pathological hallmark of Alzheimers disease (AD), and is believed to be a primary cause of dementia in elderly individuals. In the present work, we simulated the conformational evolution of A{beta}42 dimers in solution and in membrane-like environment to explore the folding of A{beta}42 along fibrillation. The molecular dynamics (MD) simulation was steered by experimental internuclear distance restraints obtained using solid-state nuclear magnetic resonance (ssNMR) spectroscopy. Our results revealed that several hydrophobic and polar motifs within the A{beta}42 sequence played key roles in the early-stage nucleation process of fibrillation and those motifs are also the stabilizing agents in the mature fibrils judged by the energy contribution. Our results also indicated that the peptide association with membrane bilayers could modulate the structural evolution pathways towards fibrillation. These findings contributed to a better understanding of the molecular level structural polymorphisms inherent to A{beta}42 fibrils. Further, the current work demonstrated that the combination of MD simulations with ssNMR-based experimental restraints provided a reliable method for studying structural changes of A{beta}. HighlightO_LIUsing solid-state NMR restraints guided molecular dynamic simulation, {beta}-amyloid dimers displayed consistent {beta}-strand-prone regions, which are major stabilizing segments for mature fibrils. C_LIO_LI{beta}-amyloid dimers evolved differently with or without interacting with the lipid bilayers. C_LIO_LIExperimental restraints guided simulation provided molecular level insights about early-stage interactions along the progress of {beta}-amyloid fibrillation C_LI

The repair of photo-induced DNA lesions through nucleotide excision repair machinery is still the source of important questions. It has been observed that the repair rate of the different cyclobutane pyrimidine dimers, i.e. the photoproducts induced by dimerization of two {pi}-stacked pyrimidines (T<>T, T<>C, C<>T, C<>C), depends on the nucleobases involved in the lesion. TT derivatives (T<>T) are removed more slowly than those containing cytosine, especially in 5. Using all-atom molecular dynamics simulations and free-energy calculations, we demonstrate that the variation of the repair rate observed in human skin and in cultured cutaneous cell is associated to the recognition of the four lesions by the DDB2 protein moiety, and more specifically by the differential structural deformation induced on the complementary strand. Indeed, while C<>C and C<>T induce a larger deviation on the groove parameters, T<>T and T<>C, instead, affect DNA structure to a lesser extent. less affected. These effects then hamper differentially the downstream recruitment of the repair complexes. The observed DNA deformation correlates with the experimental repair rate and provides a structural rationale for the different repair rates of CPD by nucleotide excision repair machinery. GRAPHICAL ABSTRACT O_FIG O_LINKSMALLFIG WIDTH=200 HEIGHT=105 SRC="FIGDIR/small/724087v1_ufig1.gif" ALT="Figure 1"> View larger version (43K): org.highwire.dtl.DTLVardef@cf6b6dorg.highwire.dtl.DTLVardef@195e35forg.highwire.dtl.DTLVardef@1829296org.highwire.dtl.DTLVardef@165baba_HPS_FORMAT_FIGEXP M_FIG C_FIG

The NAD+ dependent deacetylase sirtuin-1 (SIRT1) is known to elicit cellular defenses against aging, cancer, and other aberrant pathologies. Previous studies have identified an intrinsically disordered region of SIRT1 comprised of N-terminal residues 1-52, herein referred to as motif A, which activates SIRT1 activity, likely through intramolecular interactions. Additionally, phosphorylation of N-terminal residues Ser27 and Ser47 has been shown to be important for regulating SIRT1 activity and stability. The lack of in vitro characterization of these effects hampers our further understanding of the role of motif A in SIRT1 regulation. In this study, we elucidate the role phosphorylation plays in motif As structure as well as its regulatory effects on SIRT1 activity against Ac-p65. We find that phosphomimetic mutation at Ser27 significantly increases the activation effect of motif A towards SIRT1. This result is supported by molecular dynamics simulations of the phosphomimetics, which reveal stabilization of different transient structures for motif A depending on whether Ser27 and Ser47 have been modified. A key finding suggested by this study is that phosphorylation of S27 appears to activate SIRT1 by causing motif A, which is intrinsically disordered in the WT, to fold into an ordered structure. This conclusion is based on both the experimental findings and simulation results. These findings contribute to our understanding of SIRT1 regulation, specifically the role played by phosphorylation within the N-terminal disordered region.

Protein translation is a highly regulated process influenced by multiple factors at the initiation, elongation, and termination stages. One notable regulatory element of the ribosome is the CAR interaction surface, a three-residue motif in the structure of the ribosome composed of C1274 and A1427 of S. cerevisiae 18S rRNA (corresponding to C1054 and A1196 in E. coli 16S rRNA) and R146 of ribosomal protein Rps3. CAR is highly conserved and positioned adjacent to the amino-acyl (A site) decoding center. It establishes hydrogen bonds with the +1 codon next in line to enter the ribosome A site, acting as an extension of the tRNA anticodon and forming base-stacking interactions with nucleotide 34 of the tRNA. However, despite CARs enzymatically strategic positioning within the ribosome, its functional relationship with the A site remains poorly characterized. Using molecular dynamics (MD) simulations, we examined the interplay between the A site and CAR site, revealing sequence-dependent modulation of H-bonding and {pi}-stacking interactions within and between the two sites. These findings highlight the interplay between the A site and CAR site, suggesting a structural and functional connection between these two regions of the ribosome that may contribute to mRNA sequence-specific tuning of translation elongation. GRAPHICAL ABSTRACT O_FIG O_LINKSMALLFIG WIDTH=200 HEIGHT=91 SRC="FIGDIR/small/714784v1_ufig1.gif" ALT="Figure 1"> View larger version (22K): org.highwire.dtl.DTLVardef@1919efaorg.highwire.dtl.DTLVardef@15c4882org.highwire.dtl.DTLVardef@19c7782org.highwire.dtl.DTLVardef@16a1246_HPS_FORMAT_FIGEXP M_FIG C_FIG

O_LIThe N-terminal signal sequences of plant cytochrome P450 enzymes are recognized as critical determinants for subcellular localization and functional diversification, yet their evolutionary drivers and mechanisms remain largely unresolved. C_LIO_LIIn this study, the evolutionary trajectories of these signals were systematically decoded through the integration of the protein language model ESM-2 with phylogenetic and selection analyses. A conserved functional fingerprint was identified. This region may serve as the essential endoplasmic reticulum targeting signal and be evolutionarily decoupled from adjacent surfaces under positive selection during lineage-specific expansions. C_LIO_LIA functional-adaptive decoupling model is proposed to explain this pattern, wherein a conserved functional core is maintained while surrounding interfaces diversify. This evolutionary architecture is interpreted as the outcome of a two-step cycle: an initial phase of positive selection driving functional innovation, followed by pervasive neutral evolution that facilitates structural exploration and potentiates future adaptations. C_LIO_LIThis work demonstrates how interpretable machine learning can be integrated with evolutionary theory to reconcile neutralist and selectionist perspectives on protein evolution. A novel framework is thus provided for understanding the layered evolution of protein modules, where structural constraint, adaptive innovation, and neutral drift operate on distinct tiers to generate functional diversity. C_LI

Intercalation of small molecules between DNA base pairs affects DNA conformation, disrupting essential cellular processes including replication, transcription, and repair. We investigated conformational changes in 18-mer DNA upon intercalation of doxorubicin, SYBR Gold and YOYO-1 using extensive MD simulations. Two main patterns for the intercalation were identified: RISE-type intercalation occurs between adjacent base pairs and extends the DNA helix with decreased twist angles, while OPEN-type intercalation proceeds through base-pair opening without significant DNA extension. Kinetic analysis revealed that association rates for intercalation followed the order: first YO-moiety (mono-intercalation) > SYBR Gold > doxorubicin > YOYO-1 (bis-intercalation). Free energy landscape showed that forces at DNA termini reached up to 117 pN during stretching. Notably, base pairs adjacent to intercalators were protected from strand separation, accompanied by additional helical unwinding. MM-PBSA/GBSA analysis revealed that the driving force for intercalation is the stacking energy, and the binding affinity was highest for minor groove binding. Persistence length decreased with single molecule binding but recovered with two molecules due to their electrostatic repulsion. Mechanical properties of intercalated DNA showed position-dependence, demonstrating that multiple intercalation modes coexist in solution. The heterogeneous nature of intercalation explains why experimental measurements reflect ensemble averages rather than single binding configurations.

DNA topology is a key regulator of chromatin structure and transcription, yet its direct role in transcription factor recognition remains unclear. Here, we investigate how distinct DNA topological states modulate binding of the Saccharomyces cerevisiae bZIP transcription factor GCN4 using topologically defined plasmids. By combining, complementary biochemical approaches, including Bio-Layer Interferometry applied here for the first time to topology-dependent protein-DNA interactions, we show that DNA supercoiling directly reshapes GCN4-DNA recognition. Positively supercoiled DNA forms more stable and persistent complexes, whereas negatively supercoiled DNA retains greater conformational heterogeneity. To interpret these effects, we performed multiscale molecular simulations. Coarse-grained simulations of plasmids recapitulate the global topology-dependent trends observed experimentally, while matched minicircle models reproduce the same behaviour at the local scale. In strong agreement with experimental data, simulations reveal that DNA topology modulates the conformational ensemble of the GCN4 basic region. Overall, positively supercoiled DNA promotes a more ordered binding mode and localized protein distribution, whereas negatively supercoiled DNA supports increased structural plasticity. These findings identify DNA topology as an active determinant of transcription factor recognition and provide a multiscale framework linking global DNA mechanics to local protein-DNA interactions. Graphical Abstract O_FIG O_LINKSMALLFIG WIDTH=200 HEIGHT=113 SRC="FIGDIR/small/722604v1_ufig1.gif" ALT="Figure 1"> View larger version (51K): org.highwire.dtl.DTLVardef@18f8ba9org.highwire.dtl.DTLVardef@11a395dorg.highwire.dtl.DTLVardef@ac093borg.highwire.dtl.DTLVardef@923212_HPS_FORMAT_FIGEXP M_FIG C_FIG

Invertebrate vision relies on bistable visual pigments flipping upon photon absorption between rhodopsin and metarhodopsin states. In living butterflies, the UV-VIS absorption spectra of rhodopsin and metarhodopsin, respectively with 11-cis and all-trans isomers of 3-hydroxy-retinal (A3) chromophore, can be conveniently recorded from the eyeshine, the light reflected from the compound eye after passing twice through the light-guiding rhabdoms. * Here, a microscope coupled with a broadband LED source and a microspectrometer was used to record photorelaxations reported in eyeshine reflection spectra. Fitting temporal exponential relaxations to log-reflectance arrays yielded transient and baseline spectra that are analogous to absorbance difference and sum, respectively. Both types of spectra were subjected to singular value decomposition and to fitting of templated visual pigment absorption spectra. * The compound eye of the high brown fritillary Fabriciana adippe was exposed to a series of second-long broadband light pulses, causing photorelaxations with time constants between 40 and 120 ms that led to 80% metarhodopsin in equilibrium. The transient and baseline spectra were fitted with pigment templates, estimating the alpha peak wavelength 547-552 nm for rhodopsin and 496-501 nm for metarhodopsin. The metarhodopsin to rhodopsin alpha peak absorbance ratio 1.25-1.35 is consistent with the isosbestic wavelength at 530 nm. The second isosbestic wavelength indicates that rhodopsin beta (UV) peak absorbs more strongly than metarhodopsin below 405 nm. * Baseline spectra, which were not explicitly analysed in previous studies, enable concatenation of exposures, monitor long-term changes of pigment, and enhance the estimation of beta peak parameters. * The method can be directly used in many butterflies and could be adapted to other insects, particularly fruitflies, facilitating studies of the relation between the visual pigment spectra and the opsin sequences. Spectroscopic results can be complemented with physiologically measured photoreceptor spectral sensitivity datasets and analysed with the same global fitting procedure.

We have conducted all atom molecular dynamics simulations of POPC and DPPC lipid bilayers using AMBER Lipid21 force field with eight different water models, including SPC/E, TIP3P, TIP3P-FB, TIP4P-FB, TIP4P-Ew, TIP4P/2005, TIP4P-D, and OPC, to identify the most compatible one without any modification. A number of parameters have been computed in order to understand the structure of the lipid bilayer: Area per lipid, Isothermal compressibility modulus, average Volume per lipid, electron density profile, bilayer thickness, X-ray and neutron scattering form factors, deuterium order parameter, and radial distribution function. The estimated Area per lipid, Isothermal compressibility factor, volume per lipid and bilayer thickness are highly consistent with experimental results for the SPC/E water model, indicating its suitability with the AMBER Lipid21 force field, insted of any modification. The bilayer electron density profiles of both the lipid bilayers demonstrate a little augmentation of water penetration with respect to the membrane surface for TIP4P-D water model. However, the experimental X-ray and neutron scattering form factors are aligning well with the simulated results for all studied water models, and TIP4P-D shows better for X-ray data. The deuterium order parameter for lipid acyl chains value less than 0.25 for all observed water models, depicting their disorderness for both the lipid bilayers. The lateral diffusion and reorientation autocorrelation function of the lipid molecules in both the bilayers are computed to reveal their dynamics across all water models. In comparison to other water models, the simulated trajectories predict better structure and reasonably fair dynamic properties for the SPC/E water model. The TIP4P-Ew water model reproduces the lateral diffusion co-efficient in close agreement with experiment. Reorientational dynamics for both the lipids in the bilayers for eight different water models are observed; the presence of slow and slowest time components corresponds to the lipid axial motion (wobble motion) and Twist/Splay motions. So, in view of the overall performance of the different water models with the AMBER Lipid21 all atom force field in reproducing membrane physical properties, the SPC/E water model appears to be an optimal choice.

Single molecule methods have become prevalent tools in elucidating molecular processes across various life science fields over the past three decades, driving breakthroughs in understanding their underlying molecular mechanisms. In our study, we employed two single-molecule methods, Forster Resonance Energy Transfer (smFRET) and high-resolution optical tweezers, to investigate the transcription of Mycobacterium tuberculosis RNA polymerase (MtbRNAP) from initiation through to termination. We aim to provide a set of comprehensive biophysical tools to deepen our current understanding of MtbRNAP and its transcription factors. These experimental assays represent an important step towards unraveling the molecular dynamics and interactions that support transcription in Mycobacterium tuberculosis.

Molecular dynamics simulations provide a "computational microscope" by which molecular phenomena can be studied at atomic resolution. However, such simulations are often expensive, usually due to a combination of system size and timescale. Various enhanced sampling methods have been proposed to overcome these challenges. Despite their effectiveness, many suffer from artifacts from energetic biases guiding the simulations, or lack of effective progress coordinates. Proteins normal modes uniquely defined by their 3D fold capture their intrinsic dynamics and could provide unbiased guidance, but how to combine these modes with molecular dynamics to generate continuous, energetically unbiased pathways has been challenging. In this study, we demonstrate that conformations generated along from normal modes using adaptive anisotropic network model provide a physical, intuitive, and generalizable progress coordinate for weighted ensemble simulations, providing a boost in efficiency and a means to generate pathways for any protein system without prior knowledge.

Reactive oxygen species generated during inflammation can oxidize viral envelope lipids, with outcomes ranging from modulated infectivity to viral inactivation. For SARS-CoV-2, the molecular mechanisms by which membrane lipid oxidation influences spike protein anchoring remain poorly understood. We use all-atom molecular dynamics (MD) simulations to quantify how graded oxidation of 1-palmitoyl-2-oleoyl-sn-glycero-3-phosphocholine (POPC) affects the anchoring of the SARS-CoV-2 spike transmembrane (TM) region in an endoplasmic-reticulum-Golgi intermediate compartment (ERGIC)-like multicomponent membrane. Viral envelopes containing 0, 25, 50, 75, and 100% oxidized POPC (PoxnoPC) corresponding to 0 - 55% oxidation of all PO-type phospholipids were simulated with the spike TM helix and cytoplasmic tail embedded in a POPC/POPE/POPI/POPS/cholesterol mixture. Steered MD and umbrella sampling were used to calculate the potential of mean force (PMF) for extracting the TM+CT region along the membrane normal. Partial oxidation (25 - 75% POPC) produced reductions in the detachment barrier that were not statistically distinguishable from the native system within the sampling uncertainty, whereas full POPC oxidation lowered the anchoring free energy by about 23% (from 606 {+/-} 39 to 464 {+/-} 38 kJ mol-1), indicating that oxidation of roughly half of the glycerophospholipids can measurably weaken spike-membrane coupling. Despite this reduction, the remaining barrier (about 180kBT ) is still large, suggesting that oxidation alone may be insufficient for spontaneous spike detachment and likely acts synergistically with mechanical forces during fusion or immune engagement. Analysis of acyl-chain order parameters, area per lipid, membrane thickness, number-density profiles, and lateral lipid clustering reveals that POPC peroxidation decreases lipid order, thins and softens the bilayer, and disrupts cholesterol-stabilized clusters that refer to large cooperative lipid assemblies (>10 lipids) identified via RDF-based clustering. These oxidation-induced changes reduce hydrophobic matching around the TM helix and facilitate its extraction from the viral envelope. Our results provide a mechanistic link between lipid peroxidation, membrane nanostructure, and spike anchoring, supporting lipid oxidation for example during cold atmospheric plasma or ozone treatment as a physically grounded contributing antiviral mechanism against SARS-CoV-2.

Amyloid beta (A{beta}), one of the hallmark proteins of Alzheimers Disease (AD), aggregates into plaques that are strongly linked to cognitive decline and neuronal death. Reducing its aggregation propensity may provide a strategy to slow the progression of AD. While chirality modulation has emerged as an innovative approach to disrupt this process, research has primarily focused on alterations at the C position, often overlooking the impact of the second chiral center, such as the C{beta} atom of Threonine. Furthermore, the underlying mechanisms governing these chiral effects remain elusive. Given the intrinsically disordered nature of the A{beta} peptide, we employed temperature-replica exchange molecular dynamics (T-REMD) simulations to explore its rugged conformational landscape. We considered sequence mutations (A2T, A2V), N-terminal chirality inversion of the first six residues (A2V1-6D and WT1-6D), and alteration of the second chiral center (C{beta}) of Threonine (A2TC{beta}). By analyzing the effect size and population change induced by these mutations and chiral modulation, we concluded that the modulation at the N-termini is not confined locally but also exerts specific effects on the central hydrophobic core (CHC) region. Inspection of their free energy landscape and representative structures reveals that the protective or pathogenic effects of these variants correlate with their similarity to the wild type (WT) ensemble. Beyond these static thermodynamics analyses, a direct connection to phase transitions was made by estimating heat capacity as a function of temperature. Both analyses predict that A2TC{beta} may exert a pathogenic effect, in contrast to the protective nature of A2T. These findings offer a deeper understanding of the effects of site-specific mutations and chirality and shed light on the development of advanced therapeutic strategies for AD.

Polypeptoids (poly-N-substituted glycines) are synthetic peptidomimetic polymers with their sidechains attached to the backbone amide nitrogen rather than the -carbon in natural peptides. Peptoids display pronounced sequence-dependent conformational flexibility arising from the absence of backbone hydrogen bonding and slow cis/trans {omega}-dihedral isomerization. Despite growing interest in peptoid-based biomaterials, a coarse-grained (CG) model compatible with the modern MARTINI 3 framework is not yet available, limiting mesoscale simulation of peptoid structure and self-assembly. In this work, we develop the first MARTINI 3 compatible peptoid CG forcefield, covering 19 commonly used residue types. Extensive all-atom reference simulations employing parallel bias metadynamics (PBMetaD) were performed to ensure converged sampling of {omega}-dihedral transitions. Bonded parameters were derived from atomistic distribution functions via direct Boltzmann inversion (DBI), while nonbonded interactions were primarily adopted from the standard MARTINI 3 parameter library. The resulting CG model reproduces structural and thermodynamic properties in close agreement with all-atom simulations, while providing up to 57-fold enhanced computational efficiency. To facilitate its adoption by the research community, we have integrated all parameters and workflows to the MARTINI-based martinize2 tool, enabling automated generation of MARTINI 3 peptoid structures and topologies. This work establishes a transferable and computational efficient framework for simulating large-scale peptoid confirmations, assemblies, membrane interactions, and nanostructure formation, and supports the rational design of next-generation sequence-specific functional peptoid-based materials.